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We have explicitly shown that QCD is the color gauge invariant theory at non-zero mass gap as 
well. It has been defined as the value of the regularized full gluon self-energy at some finite point. 
The mass gap is mainly generated by the nonlinear interaction of massless gluon modes. All this 
allows one to establish the structure of the full gluon propagator in the explicit presence of the mass 
gap. In this case, the two independent general types of formal solutions for the full gluon propagator 
as a function of the regularized mass gap have been found. The nonlinear iteration solution at which 

' the gluons remain massless is explicitly present. The existence of the solution with an effective gluon 

I mass is also demonstrated. 
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QQ ■ I. INTRODUCTION 

(N : 

' Quantum Chromodynamics (QCD) [U, 0] is widely accepted as a realistic quantum field gauge theory of the strong 
, interactions not only at the fundamental (microscopic) quark-gluon level but at the hadronic (macroscopic) level as 
rS 1' well. This means that in principle it should describe the properties of experimentally observed hadrons in terms 
JL ' of experimentally never seen colored quarks and gluons (the color confinement phenomenon), i.e., to describe the 
hadronic world from first principles - an ultimate goal of any fundamental theory. 
\ The Lagrangian of QCD, however, does not contain explicitly any of the mass scale parameters which could have 
a physical meaning even after the corresponding renormalization program is performed. This clearly shows that it is 
not enough to know it in order to calculate the physical observables in low-energy QCD from first principles. It is 
also important to know the true dynamical structure of the QCD ground state especially at large distances, which 
' may be source of the above-mentioned mass scale parameter. If it will survive the renormalization program, then 
\^ ' QCD is a complete and self-consistent theory without the need to introduce some extra degrees of freedom in order 
Cn| to generate it. In this way it may become a mass gap so needed in non-perturbative (NP) QCD in order to explain 
C^l the above-mentioned color confinement and other NP effects @. It will be responsible for the NP QCD dynamics as 
Age D is responsible for the nontrivial perturbative QCD dynamics (scale violation, asymptotic freedom ( AF) [l|, &)■ 
The propagation of gluons is one of the main dynamical effects in the QCD vacuum. In our previous work [j| 
0^ , it has been shown that the only place when the mass gap may appear is the corresponding Schwinger-Dyson (SD) 
' equation of motion for the full gluon propagator. It should be complemented by the corresponding Slavnov- Taylor 
J> ] (ST) identity (see next section). The importance of this equation is due to the fact that its solutions reflect the 
quantum-dynamical structure of the QCD ground state. It is highly nonlinear (NL) equation, and therefore the 
. , number of independent solutions, which should be considered on equal footing, is not fixed a priori. The color gauge 
. structure of this equation has been investigated in detail in the above-mentioned paper [3| . We have explicitly shown 
■ " ' that the color gauge invariance of QCD is consistent with the mass gap, generated in the gluon sector of QCD. 

Our primary goal in this investigation is to find formal solutions for the full gluon propagator as a function of the 
regularized mass gap. However, for the sake of completeness and further clarity, it is instructive to describe briefly 
the main results of Ref. 4] in the subsequent section. 



> 
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II. THE COLOR GAUGE INVARIANCE OF QCD AT NON-ZERO MASS GAP 

QCD is a SUCi) color gauge invariant theory. As underlined above, its dynamical context is determined by the 
corresponding equations of motion, among which the SD equation for the full gluon propagator plays an important 
role. It can be written as follows: 
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DM = KM + Dlp[q)inp,{q- D)D,,{q), (2.1) 

where 



D%iq) = HTM + ^L^} (2.2) 

is the free ghion propagator, and ^ is the gauge-fixing parameter. Also, here and everywhere below T^^{q) = (5^,^ — 
{q^qvl<f) = S^j^u — Lf^^{q), as usual. Ilpa{q; D) is the full gluon self-energy which depends on the full gluon propagator 
due to the non-abelian character of QCD. Thus the gluon SD equation is highly NL one. Evidently, we omit the color 
group indices, since for the gluon propagator (and hence for its self-energy) they factorize, for example D'^'Kq) = 
D^i,{q)5'^^ . Diagrammatic representation of the gluon SD equation (2.1) is shown in our previous work [4], as well as 
the detail description of the full gluon self-energy Iipa{q] D). It is the sum of a few terms which are tensors, having 
the dimensions of mass squared. All these skeleton loop integrals are therefore quadratically divergent in perturbation 
theory (PT), and so they are assumed to be regularized, as discussed below. Let us note in advance that here and 
below the signature is Euclidean, since it implies — > when — s- and vice-versa. 

A. The mass gap 

Let us introduce the general mass scale parameter A^(Z)), having the dimensions of mass squared, by the subtraction 
from the full gluon self-energy its value at g = 0. Thus, one obtains 

n%{q- D) = np,(q; D) - np,(0; D) = np,(g; D) - Sp^A^iD), (2.3) 

which is nothing but the definition of the subtracted full gluon self-energy Ilp^{q; D). Contrary to QED, QCD being 
a non-abelian gauge theory can suffer from infrared (IR) singularities in the q^ ^ Q limit due to the self-interaction 
of massless gluon modes. Thus the initial subtraction at zero in the definition (2.3) may be dangerous [l|. That is 
why in all the quantities below the dependence on the finite (slightly different from zero) dimensionless subtraction 
point a is to be understood. From a technical point of view, however, it is convenient to put formally a = in all the 
derivations below, and to restore the explicit dependence on non-zero a in all the quantities only at the final stage. At 
the same time, in all the quantities where the dependence on A (which is the dimensionless ultraviolet (UV) regulating 
parameter) and a is not shown explicitly, nevertheless, it should be assumed. For example, A^(Z?) = A.^{\,a\D) 
and similarly for all other quantities. So all the expressions are regularized. For our purpose, in principle, it is not 
important how A and a have been introduced. They should be removed at the final stage only as a result of the 
renormalization program. 

By the mass gap we understand some fixed mass squared which is related to A^(I?) as follows: 

A?{D) =A^x c{D), (2.4) 

where the dimensionless constant c{D) depends on D, while the fixed mass squared does not depend on D. It 
will be called the mass gap. As the general mass scale parameter itself and constant c{D), it may depend on all 
other dimensionless parameters of the theory, namely A^ = A^(A, a, g^), where is the coupling constant squared 
and so on. In this section we will not distinguish between A^(_D) and A^, calling both the mass gap, for simplicity. 
From the subtraction (2.3) it follows that the mass gap A^, having the dimensions of mass squared, is dynamically 
generated in the QCD gluon sector. It is defined as the value of the full gluon self-energy at some finite point (see 
discussion above). It is mainly due to the nonlinear interaction of massless gluon modes. Let us remind that no 
truncations/approximations/assumptions/, as well as no special gauge choice are made for the regularized skeleton 
loop integrals, contributing to the full gluon self-energy 

B. The transversality of the full gluon self-energy 

Let us continue with the general decompositions of the full gluon self-energy and its subtracted counterpart, which 
enter the subtraction (2.3), as follows: 
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npAq;D) = Tp,{q)q^niq^;D)+qpq,Iliq^;D), 

n;^{q;D) = Tp„{q)q^Il%q^;D) + qpq,fl%q^;D), (2.5) 

where all the invariant functions of q^ are dimensionless ones, while in addition the invariant functions 11* (g^; D) and 
^{q^-jD) cannot have the pole-type singularities in the q^ ^ Q limit, since Ilp^{0;D) — 0, by definition; otherwise 
they remain arbitrary. 

Contracting them with qp along with the subtraction (2.3), one obtains 

U{q';D)=n%q';D) + ^^, (2.6) 

and 



U{q';D)=n^q';D) + ^^. (2.7) 

q 

It is worth emphasizing that the full gluon self-energy has a massless single particle singularity due to non-zero 
mass gap A^{D), which is of the non-perturbative (NP) origin. At the same time, its subtracted counterpart cannot 
have such a singularity, as mentioned above. In other words, this means that in the explicit presence of the mass gap 
both invariant functions of the full gluon self-energy gain additional contributions due to it (of course, not only at 
some finite subtraction point q^ = ^ Q). If the mass gap is welcome in the transversal invariant function n(g^; D), 
it is not welcome in its longitudinal counterpart Ii(q^]D), since just it violates the ST identity. Let us also note in 
advance that transversality of the full gluon self-energy and its subtracted counterpart can be achieved only in the 
formal A^(Z?) = limit (for a brief discussion of all these preliminary remarks see subsections below). So in the 
general case of non-zero A~(£') only two possibilities remain. 

(i). Both are not transversal and then 



qpTlp,{q-.D) = q„q'n{q'-D)^q„[q'W{q'-D) + /^\D)]^Q, 



qpW(q;D) - q^q^W {q^ ; D) ^ q,[q^n{q^ ; D) - /\^D)]. (2.8) 



The last inequality in the first of the relations (2.8) follows from the fact that 11'' (q^; D) cannot have a single particle 
singularity — A^(_D)/g^ in order to cancel A^(I?). 

(ii). Transversality of the subtracted gluon self-energy is maintained, i.e., W{q^; D) — and then 



qpnp,{q;D)=q,q^Il{q^;D)=q,A^{D)^0, g,n^„(g; i^) = 0, (2.9) 

Contrary to the first case, now we know how precisely the transversality of the full gluon self-energy is violated. So 
it is always violated at non-zero mass scale parameter A^(D). In this connection one thing should be made perfectly 
clear. It is the initial subtraction (2.3) which leaves the subtracted gluon-self energy logarithmical divergent only, and 
hence the invariant function 11* (g^; D) is free of the quadratic divergences, but a logarithmic ones can be still present 
in it, at any D. Since the transversality condition for the full gluon self-energy is violated in these relations, that is 
why we cannot disregard A^(Z)) from the very beginning (compare with the pure quark case considered in our initial 
work [1]). 



C. The ST identity for the full gluon propagator 

In order to calculate the physical observables in QCD from first principles, we need the full gluon propagator 
rather than the full gluon self-energy. The basic relation to which the full gluon propagator should satisfy is the 
corresponding ST identity 



q,,quDf,^{q) = iS,. 



(2.10) 
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It is a consequence of the color gauge invarince/symmetry of QCD, and therefore "is an exact constraint on any 
solution to QCD" This is true for any other ST identities. Being a result of this exact symmetry, it is the general 
one, and it is important for the renormalization of QCD. If some equation, relation or the regularization scheme, etc. 
do not satisfy it automatically, i.e., without any additional conditions, then they should be modified and not this 
identity (identity is an equality, where both sides are the same, i.e., there is no room for additional conditions). In 
other words, all the relations, equations, regularization schemes, etc. should be adjusted to it and not vice versa. It 
implies that the general tensor decomposition of the full gluon propagator is 



D^,{q) = I {T^,{q)d{q^) + ^L^,{q)} ^, (2.11) 

where the invariant function d{q'^) is the corresponding Lorentz structure of the full gluon propagator (sometimes 
we will call it as the full effective charge ("running"), for simplicity). Let us emphasize once more that these basic 
relations are to be satisfied in any case, for example, whether the mass gap or any other mass scale parameter is put 
formally zero or not. 

On account of the exact relations (2.5), (2.6) and (2.7), the initial gluon SD equation (2.1) can be equivalently 
re- written down as follows: 



D^.iq) - DIM + DlM^TpAl)[q^'^'{q^\ D) + ^^{D)]D,,{q) + DI^^L pAq)q^tl{q^ ■ D)D,,{q)- (2-12) 

Contracting this equation with and qu^ one arrives at qi_iqi, D^Aq) =i^- iS,^n{q^]D), so the ST identity (2.10) 
is not automatically satisfied. In order to get from this relation the ST identity, one needs to put Il{q^;D) — 0, 
which is equivalent to Il''{q^;D) — — (A^(I?)/q^), as it follows from the relation (2.6). This, however, is impossible 
since n*(q^;D) cannot have the power- type singularities at small q^, as underlined above. The only solution to the 
previous relation is to disregard A^(D) from the very beginning, i.e., put formally zero A'^{D) = everywhere. In 
this case from all the relations it follows that the gluon full self-energy coincides with its subtracted counterpart, and 
both quantities become purely transversal, i.e., Il{q^;D) = W{q^;D) and Yl{q^;D) — Il''{q^;D) = (see relations 
(2.5)-(2.7)). 

The one way to satisfy the ST identity and thus to maintain the color gauge structure of QCD is to discard the 
mass gap A^(_D) from the very beginning, i.e., put it formally zero A^(D) = in all the equations, relations, etc. In 
this limit the initial gluon SD equation (2.12) is modified to 

D^Jiq) = D%iq) + DlM^T,M)q^IV^{q^- DP^)D^Jiq), (2.13) 
and the corresponding Lorentz structure which appears in Eq. (2.11) becomes 



It is easy to see that the gluon SD equation (2.13) automatically satisfies the ST identity (2.10) now. Evidently, in the 
formal A?{D) — limit we denote D^^{q) and d{q^) as D^J {q) and d^'^{q^), respectively (for reason see below). As 
it has been pointed out in Ref. in this case there will be no problems for ghosts to accomplish their role, namely 
to cancel the longitudinal component in the full gluon propagator (2.13). 



D. The general structure of the full gluon propagator 

The formal A^(Z?) — limit is a real way how to preserve the color gauge invariance in QCD. Then a natural 
question arises why does the mass gap A^(£)) exist in this theory at all? There is no doubt that the color gauge 
invariance of QCD should be maintained at non-zero mass gap as well, since it is explicitly present in the full gluon 
self-energy, and hence in the full gluon propagator. However, by keeping it "alive" , the two important problems arise. 
The first problem is how to replace the original gluon SD equation (2.12), since it is not consistent with the ST 
identity unless the mass gap is discarded from the very beginning (see above). The second problem is how to make 
the full gluon propagator purely transversal when the mass gap is explicitly present. 

By introducing the spurious technics we were able to show that the ST identity (2.10) can be automatically satisfied 
at non-zero mass gap A^(Z?) as well. In other words, our aim is to save the mass gap in the transversal invariant 
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function (2.7), while removing it from the longitudinal invariant function (2.6), but without going formally to the PT 
= limit. In order to keep the mass gap "alive", and, at the same time, to satisfy the ST identity (2.10), we 
introduced a temporary dependence on A^(il') in the free gluon propagator, thus making it an auxiliary (spurious) 
free gluon propagator. Substituting it into the initial gluon SD equation (2.12) and restoring again the dependence on 
the free gluon propagator, such obtained gluon SD equation should satisfy the ST identity (2.10). After doing some 
tedious algebra, one finally obtains [41 



DA^) = DlM)+Dlp{q)iTp,{q)[q^W{q^-D) + ^''mDa.{q). (2.15) 

Such modified gluon SD equation (2.15) is satisfied by the same expression for the Lorentz structure d{q^) in Eq. (2.11) 
as the original gluon SD equation (2.12), namely 



which is not surprising, since the original gluon SD equation (2.12) and its modified version (2.15) differ from each 
other only by the longitudinal (unphysical) part. 

However, the important observation is that now it is not required to put the mass gap A^(Z?) formally 
zero everywhere. The spurious mechanism does not affect the dynamical context of the original gluon SD equation. 
In other words, it makes it possible to retain the mass gap in the transversal part of the gluon SD equation, and, at 
the same time, to cancel the term in its longitudinal part, which violates the ST identity. In this way, the modified 
gluon SD equation (2.15) satisfies automatically the ST identity (2.10). 

Due to AF in QCD the PT regime is realized at — > oo. In this limit all the Green's functions arc possible to 
approximate by their free PT counterparts (up to the corresponding PT logarithms). However, from the relation 
(2.16) it follows that in this limit the mass gap term contribution A^(L))/g^ is only next-to-next-to-leading order 
one. The leading order contribution is the subtracted gluon self-energy 11* (g^;!)), which behaves like Ing^ in this 
limit, as mentioned above. The constant 1 is the next-to-leading order term in the ^ oo limit. Such a special 
structure of the relation (2.16), namely the mass gap enters it through the combination A^{D)/q'^ in its denominator 
only, explains immediately why the mass gap A^{D) is not important in PT. From this structure it follows that the 
PT regime at ^ oo is effectively equivalent to the formal A^(iD) = limit and vice versa. That is the reason 
why this limit can be called the PT limit. And that is why we denote D^^{q\ A^ = 0) = D^^{q\ 0) = Dj^J{q), and 
hence d(q^; A^ = 0) = d{q^;0) = d^^{q^), etc., in accordance with the previous notations. Let us note, however, that 
sometimes it is useful to distinguish between the asymptotic suppression of the mass gap contribution A^/g^ in the 
q^ ^ oo limit and the formal PT A^ = limit (see our subsequent paper). 

Thus the formal PT A^(Z?) — limit exists, and it is a regular one. As it follows from above, in this limit one 
recovers the PT QCD system of equations (2.13)-(2.14) from the NP QCD one (2.15)-(2.16). So, we distinguish 
between the PT and NP phases in QCD by the explicit presence of the mass gap. Its aim is to be responsible for the 
NP QCD dynamics, since it dominates at ^ in the "solution" (2.16). When it is put formally zero, then the PT 
phase survives only. Evidently, when such a scale is explicitly present then the QCD coupling constant plays no role 
in the NP QCD dynamics. 



E. Transversality of the relevant full gluon propagator 



The NP QCD system of equations(2.15)-(2.16) depends explicitly on the mass gap A^. As it has been discussed in 
detail in our previous work then the ghosts are not able to cancel the longitudinal component in the full gluon 
propagator, i.e., they are of no use in this case (the transversality condition for the full gluon self-energy is always 
violated, see relations (2.8) and (2.9)). This is the price we have paid to keep the mass gap "alive" in the full gluon 
propagator. Our aim here is to formulate a method which allows one to make the gluon propagator, relevant for NP 
QCD, purely transversal in a gauge invariant way, even if the mass gap is explicitly present. 

For this purpose let us define the truly NP (TNP) part of the full gluon propagator as follows: 



Dl!:^P{q- A') = D,,,{q; A^) - D,,,{q; A^ = 0) = D^,{q; A^) - i^^J(g), (2.17) 

i.e., the subtraction is made with respect to the mass gap A^, and therefore the separation between these two terms 
is exact. So it becomes 
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j~)TNP 



iq;A^) = zT^,{q) dig'; A' 



jPT 



r 



1 



(2.18) 



where the exphcit expression for the TNP Lorentz structure d'^ ^ ^ {q^ ; A^) — d{q^;A^) 
from the relations (2.16) and (2.14) for d{q^; A^) and d^^ [q^)^ respectively. 
The subtraction (2.17) is equivalent to 



q 

d^^{q^) can be obtained 



D,Aq; A') = DiriT, A^) + D^J{q). (2.19) 

The TNP gluon propagator (2.18) does not survive in the formal FT A'^ = limit. This means that it is free of the PT 
contributions, by construction. The full gluon propagator in this limit is reduced to its PT counterpart. This means 
that the full gluon propagator, being also NP, nevertheless, is "contaminated" by them. The TNP gluon propagator is 
purely transversal in a gauge invariant way (no special (Landau) gauge choice by hand), while its full counterpart has a 
longitudinal component as well. There is no doubt that the true NP dynamics of the full gluon propagator is completely 
contained in its TNP part, since the subtraction (2.19) is nothing but adding zero to the full gluon propagator. 
We can write D^,iq;A^) = i {T^,{q)diq'; A^) + ^L^^iq)} {l/q^) - iT^,{q)dP^ {q^){l/q^) + iT^,{q)dPT {q^){l/q^) = 
Dj^^^{q; A'^) + Dj^J {q), and so the true NP dynamics in the full gluon propagator is not affected, but contrary exactly 
separated from its PT dynamics, indeed. In other words, the TNP gluon propagator is the full gluon propagator but 
free of its PT "tail". 

Taking this important observation into account, we propose instead of the full gluon propagator to use its TNP 
counterpart (2.18) as the relevant gluon propagator for NP QCD, i.e., to replace 



D,^{q; A') ^ Dl^P{q- A') = D,^{q; A') - D^^iq), (2.20) 

and hence d{q^; A^) ^ d^^^(g2; A^) = ^(g^; A^) - d^'^iq^). 

The subtraction (2.20) plays effectively the role of ghosts in our proposal. However, the ghosts cancel only the 
longitudinal component in the PT gluon propagator, while our proposal leads to the cancellation of the PT contri- 
bution in the full gluon propagator as well (and thus to an automatical cancellation of its longitudinal component). 
Nevertheless, this is not a problem, since the mass gap is not survived in the formal PT limit, anyway. 

In fact, our proposal is reduced to a rather simple prescription. If one knows a full gluon propagator, and is able 
to identify the mass scale parameter responsible for the NP dynamics in it, then the full gluon propagator should 
be replaced in accordance with the subtraction (2.20). The only problem with it is that, being exact, it may not be 
unique. However, the uniqueness of such kind of separation can be achieved only in the explicit solution for the full 
gluon propagator as a function of the mass gap (see below). Anyway, this subtraction is a first necessary step, which 
guarantees transversality of the TNP gluon propagator Dj^^^ {q; A^) without losing even one bit of information on 
the true NP dynamics in the full gluon propagator D|^^^{q^, A^). At the same time, its non-trivial PT dynamics is 
completely saved in its PT part D^^J^q). So it is worth emphasizing that the both terms in the subtraction (2.19) 
are valid in the whole momentum range, i.e., they are not asymptotics. 

The full gluon propagator (2.19), keeping the mass gap "alive" , is not "physical" in the sense that it cannot be made 
transversal by ghosts. Therefore it cannot be used for numerical calculations of the physical observables from first 
principles. However, our proposal makes it possible to present it as the exact sum of the two "physical" propagators. 
The TNP gluon propagator is automatically transversal, by construction. It fully contains all the information of the 
full gluon propagator on its NP context. Just it should be used in accordance with the prescription (2.20) in order 
to calculate the physical observables in low-energy QCD. In high-energy QCD the PT gluon propagator (2.13) is to 
be used. It is free of the mass gap and the ghosts can cancel its longitudinal component, making it thus transversal 
("physical"). 

Concluding, in this section we have briefly remind how to preserve the color gauge invariance/symmetry in QCD 
at non-zero mass gap. This means that from now on we can forget the relations (2.8) and (2.9) at all, since there 
are no any more their negative consequences for the truly NP QCD. In this connection let us remind the initial 
subtraction (2.3) has been done in a gauge invariant way (i.e., not in a separate propagators, which enter the skeleton 
loop integrals, contributing to the full gluon self-energy). 



III. MASSIVE SOLUTION 



One of the direct consequences of the explicit presence of the mass gap in the full gluon propagator is that the 
gluon may acquire an effective mass, indeed From Eq. (2.16) it follows that 
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where instead of the dependence on D the dependence on ^ is exphcitly shown, while here and below the dependence 
on all other parameters is not shown, for simplicity. The full gluon propagator (2.11) may have a pole-type solution 
at the finite point if and only if the denominator in Eq. (3.1) has a zero at this point = —mi (Euclidean signature), 



-ml-mln%-ml;0 + A^c{O^Q, (3.2) 

where = m'^(^) is an effective gluon mass. The previous equation is a transcendental equation for its determination. 
Evidently, the number of its solutions is not fixed, a priori. Excluding the mass gap, one obtains that the denominator 
in the full gluon propagator becomes 

+ q'n%q'; + A'c{0 = + + g^ff ^) + m^n^ (^m^- C). (3.3) 



Let us now expand W {q^ ; ^) in a Taylor series near 



™9- 



Wiq^;0 = n^(-K; + {<f + mlWi-ml- + 0[{q^ + m^f) . (3.4) 
Substituting this expansion into the previous relation and after doing some tedious algebra, one obtains 

q^ + m] + q^ff (g2; + (-™^; = {q^ + ml)[l + (-m^; ^ - m^n'^(-m2. ^ ■^(q^; ^)], (3.5) 

where n'*'^(q'^; ^) = at = ~'^o ^-nd it is regular at small (f'] otherwise it remains arbitrary. 
The full gluon propagator (2.11) thus now looks 



Z^{ml) . 1 

'((z2+m2)[l + ff^«(g2;m2)] ^ ^^^a^''^'^)^ 



D^,:.{q;ml) = iT^^iq)-^ J :r-r + i^Lf,^{q)^, (3.6) 



where, for future purpose, in the invariant function 11^'^ {q^;m^) instead of the gauge-fixing parameter ^ we intro- 
duced the dependence on the gluon effective mass squared m^, which depends on ^ itself. The gluon propagator's 
renormalization constant is 

= TTTiTT 2~EV~ 2?i7I7 2~E\- 

" 1 -I- 1P(— m^; 4) — m^li'*(— rn^; 4) 

In the formal PT limit = 0, an effective gluon mass is also zero, rn'^{£^) = 0, as it follows from Eq. (3.2). So an 
effective gluon mass is the NP effect. At the same time, it cannot be interpreted as the "physical" gluon mass, since it 
remains explicitly gauge-dependent quantity (at least at this stage). In other words, we were unable to renormalize it 
along with the gluon propagator (3.6). In the formal PT = mg(^) — limit the gluon propagator's renormalization 
constant (3.7) becomes the standard one 0,13], namely 

It is interesting to note that Eq. (3.2) has a second solution in the formal PT A^ = limit. In this case an effective 
gluon mass remains finite, but 1 + n'*(— to^;^) =0. So a scale responsible for the NP dynamics is not determined by 
an effective gluon mass itself, but by this condition. Its interpretation from the physical point of view is not clear. 
The massive solution (3.6) is difficult to use for the solution of the color confinement problem, since it is smooth 
in the q^ limit. However, its existence shows the general possibility for a vector particles to acquire masses 
dynamically, i.e., without so-called Higgs mechanism 6], which requires the existence of not yet discovered Higgs 
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particle. Apparently, it can be also useful in the generalization of QCD to non-zero temperature and density [7|, |8| 
(and references therein), when the gluons may indeed acquire effective masses. The above-mentioned possibility is 
due only to the internal dynamics and symmetries of the corresponding gauge theory. 

The general procedure described above in subsection E of section II can be directly applied to the massive solution 
(3.6). So it becomes 



D^M- ml) = Dl^^iq; ml) + D^Jiq), (3.9) 

where 



D 



TNPi 



(g2 



Uml) 



])[l + W'R{q^ 



q^[\ + W'R{q^:M 



(3.10) 



and 



Tt,y{q)-. 



[l + ff^«(g2;0)] 



^Lf,y{q) 



q- 



(3.11) 



Let us remind that in the massive solution the role of the mass gap is played by an effective gluon mass, so the formal 
PT limit is mi = 0. In accordance with our prescription (2.20), we should finally replace the full gluon propagator 



(3.6) as follows: D^i,{q]m ) 



r)TNP 



{q;ml), where the latter is exphcitly given in Eq. (3.10). 



IV. GENERAL NL ITERATION SOLUTION 



In order to find another type of the general formal solution for the full gluon propagator (2.15), let us begin again 
with its "solution" (2.16) which is 



diq^) . diq^; A^) ^ i^n-(^.,)^,(,)(^.;^.) . (4.1) 

where the dependence on D is replaced by the equivalent dependence on d and the relation (2.4) is already used. It 
is worth reminding that the invariant function 11" (q^; d) and c{d) are, in fact, the sum of the corresponding skeleton 
loop integrals (see section II and our initial paper [J|). Let us introduce further the dimensionless variable z = A^/g^. 
The full Lorentz structure (4.1) regularly depends on the mass gap, and hence on z. Thus it can be expand in a 
Taylor series in powers of z around zero z as follows: 



d(g2;A2) = d(g2;z) = ^zVfc(g'), (4.2) 

fc=0 

where the functions fkiq^) are the corresponding derivatives of d{q^\z) with respect to z at z = 0, which is equiv- 
alent to the PT A2 = limit. For example, /o(g2) ^ d{q^;z = 0) = d^'^ [q^) = [1 + Wiq^-d^"^]-^, fi{q^) = 

iddiq^;z)/dz),_^o = [d[l + W{q^;d) + cid)zr^/dz]^^^ = - [l + W{q^;dPT)]-\idP^) = -[d^^ iq^)]^c{dPT), and 
so on, i.e., fk{q^) — {—l)^d^'^{q^)[d^'-''{q^)c{d^'^)]^. Fortunately, these explicit expressions play no any role in what 
follows. In any case, they depend on the unknown, in general, quantities W{q^;d) and c(d), which by themselves NL 
depend on d and finally on and c(d^^). So our expansion (4.2) is nothing but the NL iteration series in powers of 
the mass gap (for the direct NL iteration procedure with d'^^ = 1 as input information see appendix A). To use also 
unknown functions fkiq'^) much more convenient from the technical point of view. However, it is worth emphasizing 
that, contrary to the relation (4.1), the expansion (4.2) can be considered now as a formal solution for d{q^), since 
fk{q^) depend on d^'^{q^), which is assumed to be "known". 

The functions fk{q^) are regular functions of the variable g^, since they finally depend on d^'^{q^) which is a regular 
function of q^ . Therefore they can be expand in a Taylor series near = (here we can put the subtraction point 
a = 0, for simplicity, since all the quantities are already regularized, i.e., they depend on a and so on, see appendix A). 
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Introducing the dimcnsionlcss variable x = q^/A/P, whore is some fixed auxiliary mass squared, it is convenient 
to present this expansion as a sum of the two terms, namely 



k 

fk{q') = Yl ^"An(0) + x'+'Bk{x), (4.3) 

n—O 

where the coefficient fkn{0) are the corresponding derivatives of the functions /fe(q^) = fkix) with respect to x at 
X = 0. Of course, these coefficients depend on the parameters of the theory such as A, a,^, (/^, and so on, which are 
not shown explicitly. The dependence on these parameters will be restored at the final stage of our derivations. The 
dimensionless functions B^.{x) are regular functions of x; otherwise they remain arbitrary. 
So the general Lorentz structure (4.2) becomes 



oo oc k 

d{q') = Yl ^'^(2;) = E ( E ^'Vfcn(O) + x"+'Bk{x)) . (4.4) 

fe=0 fe=0 n=0 

Omitting all the intermediate tedious derivations (which, nevertheless, are quite obvious), these double sums can be 
equivalently present as the sum of the three independent terms as follows: 



d{q^) = zY^'Y'^km{0) + aY<^'Yl ^>^rn{x) + d""^ {q% (4.5) 

where the constant a = xz = A^/M^ and the dimensionless functions Akm{x) are regular functions of x: otherwise 
they remain arbitrary. d^^{q^) denotes the terms which do not depend on the mass gap at all, i.e., it is nothing 
but the Lorentz structure of the PT gluon propagator (2.14), indeed. The summation over m explicitly shows that all 
iterations invoke each NP IR singularity labeled by k in the first term of the expansion (4.5). Thus it is the general 
NL formal expansion in powers of the mass gap (this is explicitly seen from appendix A). 
Going back to the gluon momentum variable q^, one obtains 



d(g2; A^) = d^^^(g^ A^) + d^^(g^) = ci^^^(g2; A^) + ci^^^(g^ A^) + ^^^(g^), (4.6) 

where the superscripts "INP" and "MPT" stand for the intrinsically NP and mixed PT parts of the TNP term, 
respectively (for reasons see discussion below). In other words, in the general NL iteration solution the TNP part 
itself is a sum of the two independent terms, i.e., (F^^{q'^;A'^) = d^^^{q^\A'^) + d^^"^ {q'^; A"^). Their explicit 
expressions are 



and 



^"'^(«';^') = (^) n^)M,') = (^)f (^)' £ ^-('")- (") 

k=0 k=0 m=0 

Here and everywhere below all the quantities depend on the parameters of the theory, namely A^ = A'^{X,a,^,g^) 
and Ak{q^) = ^^^=0 -^kmio^'i \ oi, ^, g^). At the same time, '^km depends in addition on the parameter a as well, i.e., 
*fem = ^km{\0L,l,g'^,a). 



A. The exact structure of the NL iteration solution 



The full gluon propagator (2.11) thus becomes the sum of the three independent terms, namely 
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D,Aq; A^) = DiriT, A^) + D^Jiq) = Dj^^iq; A^) + Df^i^^{q; A^) + D^J{q), (4.9) 

where 



Dinq;A^) = zr^.(g)d^^^(g2; A2)l = ^r,.(5)^L(g2; A^) (4.IO) 

with 



^('^^;A^)-E(^) *^--E(^) E'J>-' (4-11) 

fe=0 ^ fc=0 ^ m=0 

while 



A^) = zT,.(<z)d*'^^^(g2; A^)! (4.12) 



with d^'^^ (g , A^) given in Eq. (4.13) and 



(4.13) 

with d^'^{q^) given in Eq. (2.14). For the direct NL iteration procedure see appendix A, as mentioned above. 

Let us emphasize that the general problem of convergence of formal (but regularized) series, which appear in these 
relations, is irrelevant here. In other words, it does not make any sense to discuss the convergence of such kind of 
series before the renormalization program is performed (which will allow one to see whether or not the mass gap 
survives it at all). The problem how to remove the UV overlapping divergences Q and usual overall ones [l, 2, 10, li| 
is a standard one, i.e., it is not our problem, anyway (let us remind that the mass gap does not survive in the PT 
g^ — > 00 limit). Our problem will be how to deal with severe infrared (IR) (g^ — > 0) singularities due to their novelty 
and genuine (intrinsic) NP character (in this limit the mass gap dominates the structure of the full gluon propagator) . 
Fortunately, there already exists a well-elaborated mathematical formalism for this purpose, namely the distribution 
theory (D T) Il2ll. into which the dimensional regularization method (DRM) 13] should be correctly implemented (see 
also Refs. |14l.fl5l|). 

The INP part of the full gluon propagator is characterized by the presence of severe power-type (or, equivalently, 
NP) IR singularities [q^y^"^ , /c = 0, 1, 2, 3, .... So these IR singularities are defined as more singular than the power- 
type IR singularity of the free gluon propagator (g^)^^, which thus can be defined as the PT IR singularity. The 
INP part of the full gluon propagator (4.10), apart from the structure (A^/g'*), is nothing but the corresponding 
Laurent expansion (explicitly shown in Eq. (4.11)) in integer powers of q^ accompanied by the corresponding powers 
of the mass gap squared and multiplied by the q^-independent factors, the so-called residues $fc(A, a, ^, g^, a) — 
Sm=o '^km{^, Oi, 9^ i o) ■ The sum over m indicates that an infinite number of iterations (all iterations) of the above- 
mentioned corresponding regularized skeleton loop integrals invokes each severe IR singularity labeled by fc. It is 
worth emphasizing that the Laurent expansion (4.11) cannot be summed up into the some known function, since its 
residues are, in general, arbitrary. However, this arbitrariness is not a problem. The functional dependence, which 
has been established exactly, is all that matters (this will be explicitly shown in the subsequent paper). Let us note 
that the expansions (4.10)-(4.11) have been independently obtained in Ref. [T3| in a rather different way. 

The MPT part of the full gluon propagator (4.12), which has the power-type PT IR singularity only, remains 
undetermined, but depends on the mass gap (that is why we call this term as the mixed PT contribution, but it 
vanishes in the formal PT A^ = limit). This is the price we have paid to fix exactly the functional dependence of 
the INP part of the full gluon propagator. With respect to the character of the IR singularity it should be combined 
with the PT gluon propagator, leading to the so-called general PT (GPT) term, namely 



T^.(<z)d- (g^)+^L^.(g) 



PT I 



D%^^{q- A') ^ Dt'nq; A^) + D^J (q) = 



T^.(g)d^^^(g^A^) + eV(g) 



(4.14) 
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where dP^'^ {q^; = d*^^^(q^;A2) + d^^ {q^) is regular at small q^, while d^^ {q^ ; A'^ = 0) = and hence 
dGPT^q2.^2 = 0) = d^^(q^). Thus both terms MPT and PT present the PT-type contributions to the full gluon 
propagator (4.6). It is worth reminding that all the three terms, which appear in the right-hand-side of Eq. (4.9) 
are valid in the whole energy/momentum range, i.e., they are not asymptotics. At the same time, we have achieved 
the separation between the terms responsible for the NP (dominating in the IR (g^ — > 0)) and the nontrivial PT 
(dominating in the UV (g^ oo)) dynamics in the true QCD vacuum. The structure of this solution shows clearly 
that the deep IR region interesting for confinement and other NP effects is dominated by the mass gap. In the formal 
PT A'^ = limit, the nontrivial PT dynamics is all that matters. 



V. INP GLUON PROPAGATOR 



In accordance with our prescription, one should subtract all the types of the PT contributions in order to get the 
relevant gluon propagator for the truly NP QCD. As it follows from discussion above, in the case of the NL iteration 
solution, we should subtract the two terms. Doing so in Eq. (4.9), on account of Eq. (4.14), one finally obtains 

D,M, A^) A^) = D,Aq; A^) - DfJ^{q; A% (5.1) 

and hence d{q'^) d^^^{q^) as well, so that 

Dl^^iq; A^) = ^T,Aq)^L{q'; A^) = ,r,.(g)^ ^ ( — ) (5.2) 

where A^ = A^(A,a,^,5^) and $fc = $fc(A, a, ^, 5^) = ^^^Q^km{^,ct,^, g"^). In this connection, let us note that 
after the subtraction (5.1) is completed we can put the intermediate parameter a — 1, to equate thus the auxiliary 
fixed mass to the mass gap itself, i.e., put M'^ — A^, not losing generality. In the deep IR region {q^ 0) the mass 
gap is only one that's really matters. All other masses introduced from a technical point of view in order to clarify 
the derivations play only auxiliary role. 

It is important to emphasize that the INP gluon propagator (5.2) is uniquely defined because there exists a special 
regularization expansion for severe (i.e., NP) IR singularities, while for the PT IR singularity such kind of expansion 
does not exist at all (see Refs. [iJ: \lM ^^'^ references therein). This just determines the principal difference 
between the NP and PT IR singularities. It is also exactly defined because of its two features. The first one is that 
the INP gluon propagator depends only on the transversal degrees of freedom of gauge bosons. The second one is 
that in the formal PT A^ = limit the INP gluon propagator vanishes. Thus, one can conclude that the presence 
of severe IR singularities only is the first necessary condition, while the regular dependence on the mass gap and 
transversality is only second sufficient condition for the unique and exact separation of the INP gluon propagator 
from the PT gluon propagator. At the same time, the TNP gluon propagator is not uniquely defined, since it contains 
the MPT part, see Eq. (4.9). In other words, the INP gluon propagator is free of all the types of the PT contributions 
("contaminations"). Just it should replace the full gluon propagator in order to calculate the physical observables, 
processes, etc. from first principles in low-energy QCD after the corresponding renormalization program is performed. 

The INP gluon propagator satisfies its own equation of motion. For the sake of completeness, let us begin with the 
SD equation for the TNP gluon propagator [3], namely 

DiriT,^') = Dliq)^T,,iq)[^q'Wiq^■,Dn+<l''n^(.<l'■,D) + A^]D^Jiq) 

+ Dl^{q)tT,,{q)[q'U%q';D) + A']D^riQ;^') (5-3) 

with 



(9) = * {rMd^^'^iq') + ^LM} ^. (5.4) 

Here and below we omit the dependence on the mass gap in the propagators and their Lorentz structures, for simplicity. 
On account of this decomposition, the "solution" of the previous equation is 

,TNP( 2^ ^ W{q^;Dn-W{q^;D)-iAyq^) 

^'^ ' [i + n^{q^;D) + {Ayq^)][i + n-{q^;DPT)y 
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This expression coincides with the definition of (P'^^{q^) = d{q^) — d^^{q^) on account of the exphcit expressions 
(2.14) and (2.15), as it should be. 
From Eq. (4.9) it follows that 

Dir{<i)-Diri<l)+D^/^iq), (5.6) 
and substituting it into Eq. (5.3), one obtains the SD equation for the INP gluon propagator, namely 

Dir{q)= D^nq) + Dl{q)iT,Aq)[q'U%q';D) + A']D'J/^{q) 

+ Dl^iq)iTpAq)[-q^^'iq';D''^) + q^^'iq';D) + A']D!:jiq) 

+ Dl{q)^T,,{q)[q^I^^{q'■,D) + A^]Di^^{q). (5.7) 
Using the decompositions (2.2), (4.12) and (4.13), it can be simplified to 

q'Diriq)-^ - ^T,Aq){l + U^{q'■,D) + {A'/q'))d''^^(q) 

~ ^T,^{q)(-W{q';D^^) + n%q^;D) + {A^ / q^)) d^^ (q) 



~ T,^{q)T,^{q)[q'n^{q';D) + A']Di^P{q), (5.8) 

where d^^^'^{q) and d^^ {q) are given in Eqs. (4.8) and (2.14), respectively. This equation is of no practical use 
due to its complicated structure. Fortunately, we already have the explicit expression for the INP gluon propagator 
(4.10)-(4.11) or, equivalently, (5.2). It is only one to be used in order to derive renormalized gluon propagator with 
the correct confinement properties. 

However, from Eq. (5.8) it follows an important observation that like the TNP SD equation (5.3) this equation 
cannot be reduced to the free gluon propagator, when the interaction is to be switched off (i.e., setting formally 
n*(g^; D^'^) = n*((7^; D) = A^ = 0). Evidently, to the same conclusion one comes from the explicit expressions (4.8) 
and (5.5), on account of the relation d^^^iq^; A^) = d^^^iq^iA^) - d'''^^ {q^; A^), which follows from Eq. (4.6). So 
in INP QCD the gluon propagator is always "dressed" as well, and thus this theory has no free gluon propagator in its 
formalism. As it has been argued in our initial work [3| , it makes it possible to suppress the emission and absorbtion 
of the colored dressed gluons at large distances by the renormalization of the mass gap. Both the suppression of 
the dressed gluons and the absence of the free gluons are necessary for the explanation of gluon confinement by INP 
QCD (see our next paper). On the other hand, the full gluon propagator (2.19) which satisfies Eq. (2.15) is reduced 
to the free gluon propagator when the interaction is switched off. There is no mechanism to suppress the emission 
and absorbtion of the free gluons at large distances [4]. That is why the full gluon propagator (2.19) is not confining, 
while the INP one (5.2) can be. 

The subtraction (5.1) seems to be necessary, indeed. It makes the relevant gluon propagator (5.2) 
transversal and excludes the free gluons from the theory at the same time. 



VI. CONCLUSIONS 



The structure of the full gluon propagator in the presence of the regularized mass gap has been firmly established. 
We have shown explicitly that in its presence at least two independent and different typed of formal solutions for the 
regularized full gluon propagator exist. No truncations/approximations/assumptions are made in order to show the 
existence of these general types of solutions. Also, our approach, in general, and the above-mentioned solutions, in 
particular, is gauge-invariant, since no special gauge has been chosen. Let us emphasize that before the renormalization 
program is performed the gauge invariance should be understood in this sense only. 

In the presence of the mass gap the gluons may acquire an effective gluon masses, depending on the gauge choice 
(the so-called massive solution (3.6)), but a gauge-fixing parameter remains arbitrary, i.e., a gauge is not fixed by 
hand (see remarks above). Its relation to the solution of the color confinement problem is not clear, even after the 
renormalization program is performed. 

The general NL iteration solution (4.9)-(4.13) for the full gluon propagator depends explicitly on the mass gap. It 
is always severely singular in the q^ ^ limit, so the gluons remain massless, and this does not depend on the gauge 
choice. However, we argued that only the INP gluon propagator (5.2) is to be used for the numerical calculations of 
physical observables, processes, etc. in low-energy QCD from first principles. 
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It is worth emphasizing that there exists only one general restriction on the behavior of W{q'^; D), which enters the 
corresponding gluon SD equation (2.15), in the explicit presence of the mass gap within our approach, namely 



q^U'{q^;D)-.0, q^^O, (6.1) 

at any D. It stems from the second of the exact decompositions (2.5), since the subtracted gluon self-energy in 
this limit (or more precisely at q^ fj?) should go to zero. Otherwise the invariant function W{q'^;D) remains 
arbitrary (but it is logarithmic divergent at infinity). Both general types of formal solutions the massive solution 
and the NL iteration one satisfy it. The existence of some other solution(s) for the full gluon propagator, satisfying 
the general condition (6.1), should not be excluded a priori. Let us remind that the gluon SD equation (2.15) is 
highly NL, so the number of independent solutions is not fixed. Any concrete solution obtained by lattice QCD 
or by the analytical approach based on the SD system of equations is a particular case of the general types (finite 
or singular at zero gluon momentum) of the formal solutions established here. They are subject to the different 
truncations/approximations/assumptions and the concrete gauge choice imposed on the invariant function Il^{q^; D), 
which, in general, remains arbitrary but satisfying the above-mentioned general constraint (6.1) within our approach 
(see, for ex amp le recent papers [l6j, [ITI 118|, 119|, (2^ [2l|, and references therein. Let us also point out Refs. 
|23l [2^ . [25l . |26| | as well, where the gluon propagator is finite and contains the mass scale parameter. However, it, 
apparently, cannot be interpreted as gluon effective mass). 

The INP solution (5.2) is interesting for confinement, but the two important problems remain to solve. The first 
problem is how to perform the renornialization program for the regularized mass gap = A^(A, a, ^, g^), and to see 
whether the mass gap survives it or not (it has been already discussed in our previous work [4]). The second problem 
is how to treat correctly severe IR singularities {q^)~'^~^ , A: = 0, 1, 2, 3, ... inevitably present in this solution (see a few 
brief remarks above in section IV). Both problems will be addressed and solved in our subsequent paper. 
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APPENDIX A: DIRECT NL ITERATION PROCEDURE 



In order to find a formal solution for the regularized full gluon propagator (2.11), on account of its effective charge 
(2.16), let us rewrite the latter one in the form of the corresponding transcendental (i.e., not algebraic) equation, 
namely 



d{q') = 1 - (9^; d) + ^c{d) d{q') = 1 - P{q';d)d{q') 

q^ 



(Al) 



where Eq. (2.4) has been already used, and instead of D an equivalent dependence on d is introduced. It is suitable for 
the formal NL iteration procedure. For future purposes, it is convenient to introduce short-hand notations as follows: 



c(d = rf(°)+rf(i)+rf(2) + ... + rf(™) + ...) = c,„ = c„(A,a,^,g2), 
U^q^;d^d^°'> +d^^'> +d^^'> +... + d^^"''> +...) = n^(g2). 



(A2) 



and 



A2 



0,1,2,3, 



(A3) 



Via the corresponding subscript m it is explicitly seen which iteration for the gluon form factor d is actually done in 
c{d), n*(g^;(i) and P{q^;d). Let us also point out that all the invariant functions W„i{q'^) can be expand in a formal 
Taylor series near the finite subtraction point a. If it were possible to express the full gluon form factor d{q^) in 
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terms of these qiiantities then it would be the formal solution for the full gluon propagator. In fact, this is nothing 
but the skeleton loops expansion, since the regularized skeleton loop integrals, contributing to the gluon self-energy 
as mentioned above, have to be iterated. This is the so-called general NL iteration solution. This formal expansion 
is not a PT series. The magnitude of the coupling constant squared and the dependence of the regularized skeleton 
loop integrals on it is completely arbitrary. 

It is instructive to describe the general iteration procedure in some details. Evidently, d'^"-' = 1, and this corresponds 
to the approximation of the full gluon propagator by its free counterpart. Doing the first iteration in Eq. (Al), one 
thus obtains 

= 1 _ Po(g2) + ... = 1 + rf(l)(g2^ ^ (A4) 

where obviously 

rfW(g2) = _Po(g2)_ (A5) 

Carrying out the second iteration, one gets 

d{q^) = 1 - Pi(g')[l + d^'\q'')] + ... = 1 + d^'\q^) + d^^\q^) + (A6) 

where 

d(^\q^) = -d^'\q^) - Pi(g2)[l - Po{q% (A7) 
Doing the third iteration, one further obtains 

d{q') = 1 - P2(g2)[l + d^'\q') + d^'\q')] + ... = 1 + d^'\q^) + d^'\q^) + S'\q^) + (A8) 

where 



d<^^\q^) = -d(i)(g2) - d(2)(g2) _ P2(g2)[l - _ P^^q^))]^ (A9) 

and so on for the next iterations. 

Thus up to the third iteration, one finally arrives at 



°° r A2 1 r r n r ■ 

d{q^) = J2 d^"'\<i') = 1 - [n^(9') + — C2J [1 - [K{<f) + — cij [1 - ng(g2) - —CO 

TO=0 ^ ^ 1 



(AlO) 



We restrict ourselves by the iterated gluon form factor up to the third term, since this already allows to show explicitly 
some general features of the NL iteration solution. 



1. Splitting/shifting procedure 
Doing some tedious algebra, the previous expression (AlO) can be rewritten as follows: 



d{q^) = 



1 - n^(«^) + \ii{q')ni{q') - ng(g^)nf (5^)n|(g^) + ... 

A2 



q- 
A4 



ng(g^)ciC2 + nf (g^)coC2 + n^(r)coci + ... 



C2 



-C1C2 



-C0C1C2 + 



(All) 
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This formal expansion contains three different types of terms. The first type are the terms which contain only different 
combinations of n*„((j'^) (they are not multiplied by inverse powers ofq^); the third type of terms contains only different 
combinations of (A^/g^). The second type of terms contains the so-called mixed terms, containing the first and third 
types of terms in different combinations. The two last types of terms are multiplied by the corresponding powers 
of Such structure of terms will be present in each iteration term for the full gluon form factor. However, any 

of the mixed terms can be split exactly into the first and third types of terms. For this purpose the formal Taylor 
expansions for 11^ (g^) around the finite subtraction point a should be used. Thus an exact IR structure of the full 
gluon form factor (which just is our primary goal to establish) is determined not only by the third type of terms, 
ft gains contributions from the mixed terms as well, but without changing its functional dependence (see remarks 
below). To demonstrate this in some detail, it is convenient to express the previous expansion (All) in terms of 
dimensionless variables and parameters introduced in section IV, namely 



A2 

n2 ' 



M2' 



ZX ■ 



M2' 



(A12) 



where is some fixed mass squared, and /x^ is the fixed point close to = o (to be not mixed up with the tensor 
index). Also, in the formal PT = limit a = as well, since is fixed. On account of the relations (A12), the 
expansion (All) becomes 



d{x) = [l-n^(a;)+nf(a;)n^(a;)-ng(a;)nf(a;)n^(x 



+ ... 



+ z 



n^(a;)ci + ILl{x)c2 - ng(a;)nf (a;)c2 - ng(a;)n^(a;)ci - n?(a;)n^(a;)co + 

ng(x)ciC2 + ni{x)coC2 + n2(x)coci 



C0C1C2 



(A13) 



Taking into account the above-mentioned formal Taylor expansions 



00 n 

W^{x) = - a)"nW(a) = J2 [^p„fca;'=a"-'=JnW(a) 



(A14) 



n=0 n=0 fe=0 

for example, the mixed term zll2{x)ci can be then exactly split/decomposed as follows: 



012:112(0;) = CiZ 



00 

n 



n^")(a) = zPx{a) + Po(a) + 02{x). 



(A15) 



n=0 fc=0 



Here and below the dependence on all other possible parameters is not shown, for simplicity. The dimensionless 
function 02(3^) is of the order x at small X] otherwise it remains arbitrary. The first term now is to be shifted to the 
third type of terms, while the remaining terms arc to be shifted to the first type of terms. All other mixed terms of 
similar structure should be treated absolutely in the same way. 
The mixed term z'^IiQ{x)c\C2 can be split as 



00 n 

cxC2z''ni{x) = C1C20' [l]Pnfca;'=a"-'=]n[,")(a) = z''P2{a) + zN,{a) + No{a) + Oo{x), (A16) 



n=0 k=0 



where the dimensionless function Oo{x) is of the order x at small x: otherwise it remains arbitrary. Again the first 
two terms should be shifted to the third type of terms, while the last two terms should be shifted to the first type of 
terms. 

Similarly to the formal Taylor expansion (A14), we can write 



00 00 n 

W^{x)W^,{x) = W^^,{x) = ^(x - arn(±,{a) = ^ ^p^kX^a^-^Ytlio^) 



ra=0 



n=0 fe=0 



(A17) 
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Then, for example the mixed term zIlQ{x)Ill{qx)c2 can be spHt as 



oo 



n=0 



/ , 1 

k=0 



n[,^^(a) = zMi(a) + Mo (a) + Ooi(x), 



(A18) 



where the dimensionless function Oqi{x) is of the order x at small x] otherwise it remains arbitrary. Again the first 
term should be shifted to the third type of terms, while other two terms are to be shifted to the first type of terms. 

Completing this exact splitting/shifting procedure in the expansion (A13), and restoring the explicit dependence on 
the dimensional variable and parameters (A12), one can equivalently present the initial expansion (All) as follows: 



d{q^) = [-^)Bi{\a,tg\a)+(^—) B2{\a,^,g\a) + (^—) BgCA, a, 5', a) + ... + ^3(9'; A^) + ... , (A19) 

since the coefficients of the above-used expansions depend, in general, on the same set of parameters: A, a, ^,(7^, a, 
etc. The invariant function d^(q^] A^) is dimensionless, and it is free of the power- type IR singularities; otherwise it 
remains arbitrary. In the formal PT = limit it survives, and is to be reduced to the sum of the first type of terms 
in the expansion (All). In other words, it is a sum of d^^^{q'^) and dF^ {q^) up to third order, which have been 
defined in section IV. The generalization to the next iterations is almost obvious, and one finally obtains expansions 
(4.9)- (4. 13) for the full gluon propagator. 

Concluding, let us underline that the splitting/shifting procedure does not change the structure of the NL iteration 
solution at small q^ . It only changes the coefficients at inverse powers of q^ in the corresponding expansion. In other 
words, it makes it possible to rearrange the terms in the initial expansion (All) in order to get it in the final form 
(A19). Also, in the ^ limit, it is legitimate to suppress the subtracted gluon self-energy in comparison with the 
mass gap term in the initial Eq. (Al). Nevertheless, as a result of the splitting/shifting procedure, which becomes 
almost trivial in this case, one will obtain the same expansion (A19) with only different residues, as just mentioned 
above. It is worth emphasizing that residues remain completely arbitrary (undetermined) in any case. 
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